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The energy conditions of mimetic -f(R) gravity theory, together with the Dolgov-Kawasaki insta¬ 
bility will be analysed. It will be shown that the condition for the stability of the theory against 
Dolgov-Kawasaki is equivalent to the standard f(R) gravity theory. The Brans-Dicke equivalence of 
the theory is also discussed in more details. We will then obtain the parameter space of the theory 
in some special forms of f{R), which violates the strong energy condition while satisfying the weak, 
null and dominant energy conditions. 


I. INTRODUCTION 

Instead of searching for an action which is invariant under the conformal transformations, one can define a new 
metric tensor which is conformally invariant by itself. Writing any gravity action with this new metric tensor will then 
results in a conformally invariant gravity theory. This was recently done by Chamseddin and Mukhanov [1] where 
they have defined a new metric tensor, made of a dynamical metric tensor and a scalar field as 

SW = -gnvdafpdpcpg 0113 . (1) 

The metric g^ is invariant under the Weyl transformation on g^ as 

t n (x)g^. (2) 

Writing the usual Einstein-Hilbert action with the above new metric leads to a gravitational theory with an extra 
degree of freedom which behaves like a dust [1]. So, one has a new pressureless matter fluid with geometric origin 
which can play the role of cold dark matter sector of the universe. The theory is dubbed “mimetic dark matter”. An 
immediate consequence of definition (2) is that the scalar field should satisfy the relation 

g^d^dvcj) = -1. (3) 

This means that the gradient of the scalar field should be timelike vector field everywhere. One can add the above 
equation into the action as a constraint through Lagrange multiplier. It turns out that both actions yield the same field 
equations [2]. Many works have been done in the context of mimetic dark matter theory, including the Hamiltonian 
analysis [3], cosmological perturbations [4], and other generalizations of the theory [5]. The constraint equation 
actually breaks the Lorentz invariance of the action by defining a preferred timelike direction in the space-time. Such 
theories has been considered vastly in the literature, known as “Einstein-aether” theories [6]. In these theories a 
timelike vector field, known as aether vector, is introduced to the theory through a Lagrange multiplier. Many works 
has been done in the context of Einstein-aether theory, including cosmology [7] and black hole solutions [8] . One of 
the interesting facts about the Einstein-aether theory is that it can be considered as a low energy limit of the Horava- 
Lifshitz gravity [9] if the aether vector is hypersurface orthogonal. In order to make the aether vector hypersurface 
orthogonal, one can define the aether as the gradient of a scalar field T as 


M vV^VaTV^T 

This changes the Einstein-aether theory to a gravitational theory with an extra scalar degree of freedom with a timelike 
gradient. Assuming that the gradient of the scalar field T has a unit norm, makes the Einstein-aether theory and the 
mimetic dark matter theory very similar. The main difference between these two theories is that the Einstein-aether 
theory has some self-interaction terms for the aether vector. The first attempt to add some kinetic terms (other than 
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the canonical kinetic term which is already present in the theory) to the action of mimetic dark matter theory was 
done in [10] by writing the scalar version of the Einstein-aether theory using the substitution = V^</. The theory 
can then be considered as a scalar Eintein-aether gravity [10]. Unfortunately, the resulting theory is equivalent to the 
projectable Horava-Lifshitz gravity which is known to have ghost and strong coupling problems. This problem can 
simply be solved by adding to the theory a potential term which breaks the shift symmetry on the scalar field <j> [ 10 ]. 
However, one can add some other, i.e. healthy, self-interactions to the mimetic gravity action which makes the theory 
different from the projectable Horava-Lifshitz gravity, and in the mean time keeps the shift symmetry on the scalar 
field. In [11], the authors have added galileon self-interactions to the mimetic action making sure that the theory does 
not suffer from the Ostrogradski instability. The resulting theory can then explain the late-time accelerated expansion 
even in the dust dominated universe. In [12], the authors have considered the term 7 (^)(D ())) 2 to the mimetic action 
which simultaneously has higher derivative terms and breaks the scalar shift symmetry. It is interesting to note that 
the Einstein-Hilbert action plus the above term with constant 7 in the context of mimetic theory can be written as 

S = J dtd 3 xNy/^q( 3 R + K i:j K ij + (7 - 1 )K 2 ), 


where I\ij is the extrinsic curvature and K is its trace. Also we have used qp v to show the metric on the spacial 
hypersurface. This ADM form of the action shows that the following action breaks the Lorentz invariant and produces 
the canonical kinetic term of the Horava-Lifshitz gravity [13]. One can also add some potential term for the scalar 
field to the action in order to explain the effect of dark energy in this model [14]. 

Beside adding higher order interaction terms to the action, one can make the geometry richer to obtain the self- 
accelerated expanding universe. One of the most interesting theories in this category is that of f(R) gravity. Con¬ 
structing the f(R) action with the effective metric (1) may cause the universe to accelerate [15]. In [16], the authors 
have considered the dynamical system analysis of such model and showed that the universe has a stable accelerated 
expanding mode. It will then be very interesting to investigate the parameter space of the theory from the energy 
condition point of view. The mimetic-/ (R) theory can be considered as a gravitational theory coupled to a geometrical 
dark matter field which corresponds to the scalar field <j>. The modification of geometry is then responsible for the 
late-time accelerated expansion of the universe. By rewriting the theory as a Brans-Dicke theory with oj = 0, the 
expansion of the universe will be encoded to a dynamical scalar field. In this sense the model can be considered as 
a gravitational theory coupled to two scalar degrees of freedom which are responsible for both dark sectors of our 
universe. 

Despite the fact that the f(R) gravity theories are very interesting and useful, they may suffer from instabilities. In 
2003, Dolgov and Kawasaki [17] showed that the higher derivative theories may have an instability dropping the theory 
to a strong curvature regime in a very short time. One can obtain a constraint on the form of f(R) to avoid such 
an instability [18]. It is very interesting that the Palatini-/(i?) gravity does not have such an instability [20]. In [19] 
the authors found an example of f(R) gravity which passes the Dolgov-Kawasaki instability. The Dolgov-Kawasaki 
instability has been studied in the various versions of higher derivative gravities, including theories with non-minimal 
coupling between matter and geometry [21]. In the present paper we will obtain the conditions on the form of f(R) to 
avoid the Dolgov-Kawasaki instability of the mimetic -f(R) theory. We will see that the result is the same as standard 
f(R) gravity, due to the fact that the dark matter sector can not be responsible for the Dolgov-Kawasaki instability. 

The main purpose of the present work is to consider the energy conditions of the mimetic -f(R) theory and make 
some constraints on the parameter space of the theory, in order to explain the late-time accelerated expansion of the 
universe (For recent literature on the energy conditions of extended theories of gravity see [22]). For this purpose, one 
should violate the strong energy condition. The standard baryonic matter does not have such a property. However, 
the present model as well as all modified gravity theories, can violate the strong energy condition by choosing the 
dynamics of the extra d.o.f properly. This obviously does not affect the healthy behavior of the baryonic matter 
as long as it is independently conserved. One should note that the weak, null and dominant energy conditions 
should not be violated. In this paper, we will obtain the conditions of violating the strong energy condition for three 
different/important choices of the function f(R), i.e. the power law, the exponential case and also the mixed case. The 
energy conditions can be inferred from the Raychaudhuri equation which is the temporal evolution of the expansion 
scalar [23]. In the case where the matter content of the universe can be expressed as a perfect fluid, one can write [24] 


• Weak Energy Condition (WEC): 

• Null Energy Condition (NEC): 

• Strong Energy Condition (SEC): 

• Dominant Energy Condition (DEC): 


Peff > 0 
Peff + Peff > 0 
Peff + %Peff > 0 
Peff - Peff > 0. 
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We will see that the energy-momentum tensor of the present model has the form of a perfect fluid. So, we will use 
the above relations for analysing the energy conditions of the theory. 

The present paper is organized as follows: In the next section we will obtain the field equations of the model and 
construct its effective energy momentum tensor. We will then discuss about the relations of the theory with the 
Brans-Dicke theory in section III. The Dolgov-Kawasaki instability analysis will be discussed in section IV. In section 
V we will obtain the energy conditions and consider some special cases for the function f(R). In section VII we will 
summarize the results and make some discussions. 


II. THE MODEL 

Let us consider the f(R) gravity action in the mimetic theory 

S = J d A x\J—g (g^v, <j>) (/ (R (g^, <j>)) + An) , (4) 

where £ m is the matter Lagrangian and the physical metric „ is related to the auxiliary metric g^ and the scalar 
field <f> as in equation (1) which yields the constraint equation (3). The field equations can be obtained by varying 
the action with respect to the metric and the scalar field <j> with the result 

\fg^ - f'R - g^uf + ( 2 / - Rf - 3 □/' + = 0 , (5) 


and 




dA 2/ - Rf - 3□/' + -T 


= 0 , 


( 6 ) 


where prime denotes derivative with respect to the argument, and T is the trace of energy-momentum tensor. 

By taking the covariant derivative of the metric field equation (5) one can obtain 

= 0, (7) 

where we have used the scalar field equation (6) and = 0, obtained by taking the covariant derivative of 

(3). So, the ordinary matter energy-momentum tensor is independently conserved in this theory. This suggests that 
the massive test particles in this theory move along the geodesics, and no extra forces is present here [21, 25]. 

From equation (6), one can obtain 


V=9 ( 2/ - Rf - 3D/' + -T ) = CV„ 


( 8 ) 


where C^ is an integration constant. It is clear from equation (8) that the vector density C^ is either zero or time-like. 
One can see from equations (5) and (6) that in the case of vanishing the theory reduces to the standard f(R) 
gravity. In this paper we will consider non-zero values of C M and compare the results with the standard f{R) gravity 
case = 0. 

Substituting equation (8) into (5), one can obtain 


r< _ r Rfj.iy , ( f 1 __ , 1 t-rj ^ _ |—|\ r/ I C^d v (t) ic\\ 

G/jj/— 2f + \2f 2 7 ^ + y/W v (9) 

The above form of the metric field equations mimic the standard f(R) gravity except the last term in the right hand 
side which is due to the presence of a constraint on the scalar field. One should mention that /' plays the role 
of modified gravitational constant which should be positive. In the following, we will impose this contraint for the 
analysis of energy condition. 

From equation (9) one can define the effective energy-momentum tensor as 


r r e f f _ 

fi 



g^ 


y 


g^y) f' + 2 


C„d„ct> 

/V=s’ 


which is conserved, = 0, by definition. 


( 10 ) 
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It is worth mentioning that the same procedure of [2] can be done in the present case to show that (4) is equivalent 
to the action 


S = J d A xV=g{f(R) + A + 1 ) + Cm) ■ ( 11 ) 

In this form, the scalar field <f> is minimally coupled to the metric and the constraint equation (3) is added to action 
through the Lagrange multiplier. To obtain the field equations one should vary the above action with respect to A, 
the metric tensor g^ and </>. The results are exactly the same as (3), (5) and ( 6 ). 

III. BRANS-DICKE EQUIVALENCE OF MIMETIC-/)#) GRAVITY 

As is well-known, the standard f(R) gravity is equivalent to the Brans-Dicke theory with the parameter w = 0. 
Let us consider the equivalence between mimetic-/ (R) gravity and Brans-Dicke theory in the absence of the matter 
Lagrangian. The Legendre transformation for the Lagrangian (4) can be written as 

S = J d 4 x^/-g(g, $(* R(g, </>) - V(*)), (12) 


where 


^ = f'(R(gA )), VW = ^R(gA) - f(R(g,<P)), 

where the prime denotes derivation with respect to the argument. The above action shows that the action (4) is 
also equivalent to the Brans-Dicke theory with parameter ui = 0. The only difference between the usual Brans-Dicke 

theory and the above action is that the physical metric is defined by the relation ( 1 ). So, we can call this theory as 

mimetic-Brans-Dicke theory. To obtain the field equations one should vary the action with respect to the d', (f> and 
g M „. The results are 

V'(tt) = R{g), (13a) 

+ (g^D - V M V„)tt+^„V(*) - ^('LG + 3D0 + 2V(T)) = 0, (13b) 

V“(0 a (*G + 3D0+2V(*))) = 0, (13c) 

together with equation (1). We note that G is the trace of Einstein tensor. In the above equations the covariant 
derivatives and also the Einstein tensor are constructed from the physical metric g 

If we use the action ( 11 ) for the mimetic-/(I?) gravity instead of (4) and performed the Legendre transformation, 
we would obtain 


S = f d A xV^g[^R - n*)) + A or + 1)], 


(14) 


where 


V = f'(R), V(*) = *R-f(R). 

It can be easily seen that the above action is the Brans-Dicke theory with ui = 0 which minimally couples to a 
constrained scalar field through the Lagrange multiplier. By varying the action with respect to the 4^, g^ and (j), one 
will obtain the same equations (13). The A field equation then gives (3). 


IV. THE DOLGOV-KAWASAKI INSTABILITY 

In this section we will obtain the Dolgov-Kawasaki criterion for mimetic-/(I?) gravity theory. We will use the 
standard procedure of [18, 19]. Suppose that the Ricci scalar differs from its background value Rq by R = Rq + R \, 
where R\ is a perturbation around Rq. Also assume that the scalar field is perturbed from its background value 
^>o(t) = t as 4>{t,x) = <f>o(t) + In order to be consistent with solar system tests, we also assume that the 

function f(R) is very close to the Ricci scalar 


f(R) =R + ei/;(R). 
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With these assumptions, one can expand the various derivatives of the function /(R) as 


/(R) — Rq + Ri + cipo + ^qRi + •••) (15) 

f'(R) = 1 + ci/jq + ciI’qRi + ..., (16) 

f"(R)=e^ + e^ , R 1 + ..., (17) 

f'"{R) = e^ / + e^"R 1 + ..., (18) 


where ipo = ip(Ro). We have also assumed that the metric is very close to the Minkowski so that one can write 

DR = -<9 f 2 R + V 2 R, V^RV^R = -( d t Rf + (dtR) 2 . 


One should note that we have two dynamical and one constraint equations. In order to study the perturbations of 
the Ricci scalar we have to take the trace of metric field equation. The trace of equation (5) vanishes after using 
the constraint equation (3). Expanding the constraint equation upto first order gives dt<f> i = 0 with implies that 
(f> i = <fi(x). So we need to consider only the scalar field equation. 

Multiplying the scalar field equation ( 8 ) by d^<j> and using the constraint equation (3), one can obtain 


2/ - Rf - 3□/' + 

2 V~9 


Expanding the above equation up to first order gives 


where we have defined 


ib’" ■ ib"' - 

Ri + 2d t R 0 ^ ir R 1 - 2 -y^-Vi? 0 .Vi?i - V 2 i?! + m 2 eff Ri = V eff , 
Vo Vo 


1 


L eff 


3eV’o 3 ^l 


- ^Rq - |(Mo - V^oV^Ro 




Veff — D-Ro — 


2 y > 0 , 1 D tP' 0 , < 
3e^o 3 V’o 


0^07)7 + ^V^RoV^Ro- 
3 % % 


(19) 


( 20 ) 


( 21 ) 

( 22 ) 


One can see that the effective potential of the mimetic theory is changed due to the presence of a constrained scalar 
field. The effective mass is dominated by its first term which implies a constraint i/j" > 0 for a theory to avoid 
tachyonic instability. We note that i / j'q = /"(Rq). So, the present theory has the criteria 


/"(Ro) > 0, 


to avoid Dolgov-Kawsaki instability (see also [16].) We note that the same condition is hold for /(R, T) theory. 
However, non-minimal couplings between the Ricci tensor and the energy-momentum tensor can modify the above 
condition [21]. In the following we will impose this constraint for analysing the energy conditions of the theory. 


V. CONSTRAINTS ON f(R) FROM ENERGY CONDITIONS 

Let us consider the flat FRW universe 

ds 2 = — dt 2 + a{t) 2 {dx 2 + dy 2 + dz 2 ) , (23) 

which is filled with a perfect fluid with energy-momentum tensor 

= (p + p)UpU v + pg^, (24) 

where p is the energy density and p is the pressure. In this case, from equation (3), one can obtain 

4>(t) = t + c (25) 

where c is an integration constant. The scalar field <j> always appears with derivative, hence the theory has a shift 
symmetry on </>, so we can set c = 0 without loss of generality. Equation ( 8 ) implies that only the temporal component 
of Cn be non-zero and given by 

a 3 ■ ... 

Cq = — (3p - p + 4/ - 12 f'(H + 2 H 2 ) + 18 Hf + 6 /') = const. 


( 26 ) 
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Here, H is the Hubble parameter and dot represents time derivative. The effective energy-momentum tensor becomes 
diagonal with components 


1 ( 2C \ 

Peff = 12 H 2 + 6 H+-[p-f- 6 Hf + J , 

Peff = -12 H 2 -6 H+jj(p+f + AHf + 2/^ , (27) 

Using the above relations, the energy conditions for our model can be written as 

• WEC: 12H 2 + 6 H+jr^p-f- 6 Hf + > 0, 

• SEC: —24 H 2 - 12 H + j 7 (^p + 3p+2f + 6 Hf + 6/' + ^j > 0, 

•NEC: j T (^p+p + 2f l -2Hf + ^\ > 0, 

• DEC: 24H 2 + 12H + j T (p-p-2f ~2f -10i7/' + ^j >0. 


The weak energy condition should be positive because it represents the energy content of the universe. However, in 
order to have an accelerated expanding universe, the strong energy condition should be violated. In the following we 
will focus on the values of the parameter space which breaks the strong energy condition, keeping the other three 
conditions positive and satisfy f > 0 and /" > 0 . 

In order to put a better constraint on the parameters of the theory, it is helpful to write the energy conditions in 
terms of the Hubble parameter, deceleration parameter, jerk and snap with the definitions 


la. la 1 "a 

IP a' S= lP~a' 


(28) 


The present day values of the deceleration parameter,the jerk and the Hubble parameter are qo = —0.81 and jo = 2.16 
[26] and Hq = 67.74 kms _ 1 Mpc _ [27]. Unfortunately, we do not have any experimental data for the snap parameter. 
In the standard f{R) gravity theory, one can use only the weak energy condition which does not have the snap [28]. 
However, in our model, due to the presence of constraint (26), one can obtain the snap parameter in terms of qo , Hq 
and jo- The constraint equation (26) can be written in terms of cosmological parameters as 


Co — 54( — / — 77 ^(po — 3po) + 2Ho(jo ~ Qo — 2) 2 /'" + -r-^o^o + 3jo + 5q 0 + so)f" + 77 -^ 0(90 — 1 )f' 


27 


108 


(29) 


where we have set ao = 1. The weak energy condition is independent of the snap parameter and can be written as 

J, (Vo - / + 677 q (1 - qo)f - 3677q (jo -q 0 - 2)f" + 2 C^j > 0. (30) 

In the other three energy conditions the snap parameter appears, which can be substituted by equation (29), with 
the result 

. SEC jr [2po - 2/ - 72Hg (jo - go - 2)/" + 2C 0 ) < 0, 

. NEC 3 ^ ^4 Po - 4/ - 12-ffo (90 - 1)/' - 144flo 4 (jo - 9o - 2)/" + 6C 0 ) > 0, 

. DEC 3 ^ [2po - 2/ - 2AHl(q 0 - 1)/' - 7277 0 4 (j 0 - 90 - 2)/" + 6C 0 ) > 0, 

where we have changed the inequality of SEC, signalling its violation. Note that by setting 167 tG' = 1 = c in this 
paper, the Hubble parameter has dimension T , where L represents the dimension of length. Also the energy density 
p , the constant Co and the function f{R) have dimensions L~ 2 . One can then write the above energy conditions as 
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• WEC i (^ 6 n m - g + 6(1 - q 0 )g' - 36(j 0 - go - 2 )g" + 2c 0 J > 0, 

• SEC ^ (l2D m - 2g - 72(jo - g 0 - 2)<?" + 2c 0 ) < 0, 

. NEC 3 ^ ^24D m - Ag - 12(g 0 - 1 )g' - 144(j 0 - q 0 - 2)g" + 6 c 0 ) > 0, 

. DEC 3 ^ (l2fl m -2 g- 24(g 0 - 1 )<?' - 72(j 0 - g 0 - 2)<?" + 6 c 0 ) > 0, 

where we have defined dimentionless quantities 


c 0 = 


Cb 

^o 2 ’ 


“ "m — — 


6i7 0 2 ’ 


r = 


R 

W 


g = 


Hi 


(31) 


In the above relations prime represents derivative with respect to r. For Q m we have used the value fl m o = 0.31 
which is the current density of baryonic plus dark matter [26, 27]. Let us study some special cases for the function 

HR)- 


A. The Case / = R + aR n 


Note that the constant a has dimension L 2 ^ 1 n \ Introducing /? = Hi"' 1 ^a, one obtains g(r) = r + (3r n . The 
energy conditions can be written as 

• WEC: 


—/3(0.29n 2 - 1.29n + l)(10.86) n + 2c 0 + 1.86 
1 + 0.09/3n(10.86) n 


• SEC: 


-/3(0.59n 2 + 0.59n + 2)(10.86) n + 2c 0 - 18 
1 + 0.09/3n(10.86) n 


• NEC: 


—/3(0.39n 2 - 1.06n + 1.33)(10.86) n + 2c 0 - 4.76 
1 + 0.09/3n(10.86) n 


• DEC: 


—/3(0.19n 2 - 1.53n + 0.66)(10.86) n + 2c 0 + 8.48 
1 + 0.09/3n(10.86) n 


and the conditions for positivity of the first and second derivatives of f(R) can be summarized as 

1 + 0.09/3n(10.86) ra > 0 & /3n{n - l)(10.86) n > 0. 

Let us first consider the case /3 > 0. The positivity of /" implies that n ^ (0,1). One can easily see from the 
above equations that for finite /3 and Co, the maximum values for all energy conditions in terms of n occur at minus 
infinity. By considering the WEC, NEC and DEC in the limit n —> — oo, it is obvious that the positivity of NEC 
implies the positivity of WEC and DEC. This can also be seen directly from figure (1). The positivity of NEC at 
minus infinity, implies that cq > 2.38. In figure (1) we have depicted the energy conditions in terms of n for Cq = 3 
and /3 = 10. We should note that larger values for [3 makes the minimum of the functions deeper and the range of 
allowed n smaller. Furthermore, the larger values for Co will raise all the functions. The maximum value for co can be 
obtained by demanding that the SEC becomes zero at minus infinity. One can then obtain co max = 9. In figure (2) 
we have plotted the allowed range of n and /3 for co = 3 which all of the conditions (32) together with the constraints 
/', f" > 0 are satisfied. One can see from the figure that for larger values of f3 the allowed range for n becomes smaller 







EC 



FIG. 1: Energy conditions for Co = 3 and /? = 10 in the case g(r) = r + /3r n . Note that larger values for /? makes the minimum 
of the functions deeper. All functions have a maximum at minus infinity for finite Co and /3. 



FIG. 2: The allowed range of parameters n and /3 > 0 in the case Co = 3 for g(r) = r + /3r n . The left/right figures correspond 
to the values n > 1 and n < 0 respectively. The line in the figure which is marked by an arrow shows that for larger values of 
/?, the allowed range for n becomes smaller. 


as was mentioned before. In summary, in the case of /? > 0 the allowed range of Co is 2.38 < Co < 9 which means that 
the ordinary f(R) gravity theory (which corresponds to cq = 0) can not violate the strong energy condition together 
with satisfying the others. 

In the case /3 < 0, considering the constraint f" > 0 implies that only the values of n in the range n £ (0,1) 
are allowed. There is an upper bound for Co in this case. Considering SEC, one can see that for values cq > 9 the 
expression will always be positive because the term containing /3 in the numerator is always positive. So, in order to 
violate the SEC, one should assume that Cq < 9. In figure (3) we have plotted the allowed range of parameters n and 
/? for co = —1, Co = 1 and Co = 0. For larger values of Co the allowed region gets smaller but the quantitative shape 
remains the same. 


B. The Case / = ae pR 

Introducing ij = a/H q and e = Uq/ 3, one can write g{r) = rje er . In this case the positivity of /' implies that ge > 0 
and the positivity of f" results in rje 2 > 0. So, both e and parameters should be positive. With this in hand, the 
energy conditions can be simplified to 

• WEC: —(34.92e 2 - 10.86e + 1 )g + (2c 0 + 1.86)e- 10 ' 86£ > 0, 

• SEC: — (69.84e 2 + 2)r? + (2c 0 + 3.72) e - ia86e < 0, 

• NEC: — (46.56e 2 - 7.24e + 1.33)?? + (2c 0 + 2.48) e - 10 ' 86e > 0, 

• DEC: —(23.28e 2 - 14.48e + 0.66)?? + (2c 0 + 1.24) e - 10 ' 86e > 0, 

In figure (4) we have plotted the energy conditions assuming r] = 1. For this choice one can show that there is no 
allowed region for e < 0.029 because the absolute value of the SEC becomes smaller than the NEC. The shaded region 
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FIG. 3: The allowed range of parameters n and fi in the case g(r) = r + /3r n for three cases Co = —1, Co = 0 and Co = 1. 


CO 



FIG. 4: The energy conditions in the case of g = 1 for g(r) = ge er . The shaded area represents the allowed range in the Co, e 
plane. 


represents the allowed range of parameters Co and e for which the SEC is violated and the others are satisfied. One 
can see that the Co < 0 region is bounded from below, which is also true for other values of 77 > 0. The reason is that 
the e polynomials in the Weak and Null energy conditions are always positive. Considering the minus sign in front of 
the polynomials and the positivity of 77 , one can see that Co should be greater that —0.93. 

In figure (5), We have plotted the allowed range of parameters e and 77 in three different cases Co = 0.2, co = —0.2 
and co = 0 . One can see from this figure that the mimetic -f(R) gravity gives a wider range of parameters compared 
to the f(R) case. A note on larger values of Co > 0 is in order here. In this case, the Weak, Null and Dominant energy 
conditions can be satisfied easily. But in order to violate the SEC one should choose small values of e and large values 
of 77 in such a way that the polynomial term dominates the exponential term. So, the qualitative shape of the figure 
for cq > 0 remains the same, and the only change is the movement of the shape closer to the 77 axis. 


VI. THE CASE f = R + aR n e pR 

By defining e = Hq and 77 = H^" ^a one can obtain g{r) = r + rjr n e er . The energy conditions in this case can be 
written as 

• WEC: ^10.866+2.38*. ( _ 34. 92e 2 _ 0 .3n 2 + 10.86e + 1.3n - 6.43en - l) + 2c 0 + 1.86 > 0, 

• SEC: ^10.866+2.38™ ( _ 6 9. 84e 2 _ 0 . 59 n 2 _ I2.86en + 0.59n - 2) + 2c 0 - 18 < 0, 

• NEC: ^10.866+2.38™ ( _ 46.56e 2 - 0.39n 2 + 7.24e + 1.061n - 8.57en - 1.33) + 2c 0 - 4.76 > 0, 

• DEC: ^gio.866+2.38™ ( _ 23.28e 2 - 0.2?r 2 + 14.48e + 1.53n - 4.29en - 0.67) + 2c 0 + 8.48 > 0, 
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FIG. 5: The allowed range of parameters for Co = 0.2 (the Solid line), Co = 0 (the Dashed line) and Co = —0.2 (the Dotted line) 
for g(r ) = ge Er . 


cO c0 



FIG. 6: The allowed range for Co as a function of n for g = —1, in the case g(r) = r + gr 71 e ~ 0 09nr . Each line represents the 
value of Co where the corresponding energy condition vanishes. Below the line, the energy condition is violated and above the 
line, the energy condition will be satisfied. 


Also the positivity of f and /" is reduced to 

^ e 10.86e+2.38n( e + Q 0 9n) _|_ 1 > 0 , (35) 

g (e 2 + O.Oln 2 + 0.18en - O.Oln) > 0. (36) 

where we have used r n = e lnr for simplifying the equations. 


A. special case e = — 0.09n 

Let us consider a special case where /' > 0 is automatically satisfied. Solving for f" > 

g > 0 n <0 or n > 5.26, 

g < 0 0 < n < 5.26. 

The above equations mimic the conditions in subsection V A with different range for n. 
that for 77 > 0, the n > 5.26 values imply the positivity of SEC. As in subsection V A, in order to obtain a constraint 
on the parameter cq one should consider the limit n —> — 00 . However, the behavior of the energy conditions at minus 
infinity is the same as in V A, which implies that 2.38 < Co < 9. 

For negative values of 77 , we have plotted the allowed range of cq as a function of n in figure ( 6 ). We have split the 
allowed range of n into two pieces because possible values of Co increase very fast as n reaches 5 due to the exponential 
term in the energy conditions. From the figure, one can see that the co = 0 case is possible in this case. The f(R) 
theory can only adopt 77 < 0 with 1.27 < n < 2.98. On the other hand the mimetic-/(i?) gravity extend the range of 
n to all values less than n = 5.26. 


0 one can find 

(37) 

However, one can see here 
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0.2 

0.0 

- 0.2 

lu 

-0.4 

- 0.6 


- 0.8 


FIG. 7: The allowed range of e and n in the case g(r) = r + r\r n e^ r . The right figure corresponds to r/ = 0.5 and the left figure 
corresponds to 77 = —0.5. In each figure, the gray area is the allowed region in which the first and the second derivative of g(r) 
is positive. The allowed range for C = 10, C = —10 and C = 0 (which corresponds to the f(R) gravity case) are shown in each 
cases. The solid line corresponds to the ordinary f(R) gravity. The dotted line corresponds to the C = 10 case and the Dashed 
line corresponds to the C = —10 case. 




B. general case 

In figure (7), we have plotted the allowed range of parameters for general e and n. The left/right figures correspond 
to negative/positive values of 77 respectively. Gray region corresponds to the allowed region of parameters in which 
the first and second derivatives of f(R) be positive. One can see that the allowed range for 77 < 0 is bounded. We 
have specified the allowed range of parameters for three different cases; negative and positive values of cq and the 
vanishing co- For the case of positive 77 , one can see that the cq > 0 values gives us more parameter space. On the 
other hand, the parameter space for co < 0 values is a subset of cq =0, which is a subset of co > 0 case. In summary, 
one can see that the mimetic f(R) gravity will enlarge the parameter space of the theory. The 77 < 0 case is more 
interesting, for which the allowed range for co > 0 and Co < 0 cases do not include the Co = 0 case. Again, we have a 
wider parameter range in this case compared to the ordinary f(R) gravity. 


VII. DISCUSSIONS AND FINAL REMARKS 

In this paper we have considered the f(R) gravity in the context of mimetic dark matter theory. The mimetic 
dark matter theory introduces a very simple procedure to construct a conformally invariant gravitational action by 
expressing the physical metric in terms of a scalar field and a dynamical metric. Surprisingly, this procedure results in 
a degree of freedom which behaves like a dust and can be considered as a geometrical candidate of dark matter. The 
action is equivalent to a gravitational theory with a constraint which is added by a Lagrange multiplier. In this point 
of view, the Lagrange multiplier plays the role of the energy density of dark matter and the covariant derivative of 
the scalar field plays the role of its velocity. One of the interesting fact about the mimetic theory is that the equation 
of motion of the metric tensor has a general form 

+ F a a = 0. 

Here = 0 is the metric field equation in the non-mimetic gravity theory. The presence of scalar field results in 
the addition of to the metric equation. In the Lagrange multiplier point of view, the trace of the metric equation 
is the Lagrange multiplier itself, which represents the energy density of the dark matter. 

The f(R) gravity theory is equivalent to the special case of mimetic-/(i?) theory with vanishing Co (see equation 
( 8 )). So, the mimetic -f(R) theory is more general, in the sense that it yields more solutions and provides a wider 
range of parameters with respect to the cases of non-vanishing Cq. 

The mimetic-/(.R) theory is equivalent to a gravitational theory with two scalar fields. One of them is responsible 
for the dark matter sector of the universe and makes the theory conformally invariant. The other is a scalar field 
which is added to the action via Brans-Dicke interaction with to = 0, and is responsible for the dark energy content of 
the universe and makes the universe accelerate. The Dolgov-Kawasaki criterion can be found to be equivalent to the 
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usual f{R) gravity theory. This can be explained by the fact that the mimetic theory only adds a dust like scalar field 
which can not change the behavior of the expansion radically. However, as we have seen in this paper, the existence 
of this scalar field can change the parameter space of the theory. 

The energy condition analysis of the theory shows that various forms of f(R) can be used to obtain the accelerating 
universe. In the case f(R) = R + aR n , we have shown that for positive values of a, the constant Co should always 
be in the range 2.38 < Cq < 9. This means that the ordinary f(R) gravity with positive coupling can not produce 
a healthy accelerating phase for our universe. However, with negative values of the coupling a one can have reliable 
accelerating cosmology for both signs of Co and also for Co = 0. In the case f(R) = a exp (/3R), we have shown that in 
order that gravity remains attractive and the Dolgov-Kawasaki instability disappears, both the parameters [3 and a 
should be positive. An interesting fact about this case is that the parameter Co is bounded from below with c m i n < 0. 
We have also considered a more general case f(R) = R + aR n exp(/3i?). We have seen that for both signs of the 
coupling a, the Co = 0 case can produce a healthy accelerating universe. In the case rj > 0, the Co > 0 values gives 
more parameter range. 

In summary, we have seen in this paper that the mimetic-/(i?) gravity generalizes the f(R) gravity theory which 
gives us more parameter space. This will produce an accelerating phase of the universe, keeping the Weak, Dominant 
and Null energy conditions positive, and remaining in the validity domain of Dolgov-Kawasaki criterion. 
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